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ABSTRACT 

A torus develops a state of suspended accretion against a magnetic wall 
around a rapidly rotating black hole formed in core-collapse hypernovae. It 
hereby emits about 10% of the black hole spin-energy in gravitational radia- 
tion from a finite number of multipole mass moments. We quantify the relation 
between the frequency of quadrupole gravitational radiation and the energy out- 
put in torus winds by fg^ ~ 470Hz(E^/4 x lO^^erg)^/^)^?]^©/]^)^/^ where 
M denotes the mass of the black hole. We propose that irradiates the rem- 
nant stellar envelope from within. We identify E^ with energies ~ 10^^ erg 
inferred from X-ray observations on matter injecta; and the poloidal curvature in 
the magnetic wall with the horizon opening angle in baryon poor outflows that 
power true GRB energies of -E^ ~ 3 x 10^^ erg. 

Subject headings: black hole physics — gamma-rays: bursts and theory - gravi- 
tational waves 



1. Introduction 



Stellar mass black holes surrounded by a torus of a few tenths of a solar mass may form in 
prompt core-collapse of young massive stars by conservation of mass and angular momentum, 
as in the Woosley-Paczynski-Brown scenario of hypernovae (Woosley 1993; Paczyhski 1998; 
Brown et al. 2000). These are transient systems which are interesting from the point of view 
of the first law of thermodynamics. A torus surrounded a rapidly rotating black hole develops 
a suspended accretion state for the lifetime of rapid spin of the black hole (van Putten 2001) 
- generally a secular timescale. The torus catalyzes some 10% of black hole spin-energy into 
gravitational radiation, with additional winds, thermal and MeV neutrino emissions (van 
Putten & Levinson 2002). Gravitational radiation is emitted by a finite number of mass 
multiple moments produced by a Papaloizou-Pringle instability (Papaloizou & Pringle 1984; 
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Goldreich et al. 1986; van Putten 2002); torus winds form analogously to pulsar winds, and 
MeV temperatures develop by internal heating due to shear in response to competing torques 
acting on the inner and the outer face of the torus. Similar considerations apply to black 
hole-torus systems formed in black hole-neutron star coalescence, provided that the black 
hole spins rapidly. 

The predicted gravitational radiation may be determined by upcoming gravitational 
wave experiments LIGO (Abramovici et al. 1992), VIRGO (Bradaschia et al. 1992), TAMA 
(Masaki ct al. 2002) and others, possibly in combination with any of the bar or sphere 
detectors presently being developed. Accurate prediction of the gravitational wave-frequency 
is important in estimating the expected signal to noise ratios (e.g.. Cutler & Thorne (2002)), 
or in the design of dedicated experiments, such as dual recycling interferometry (Harry et 
al. 2002). 

In this Letter, we quantify the relation between the gravitational wave-frequency and 
the associated energetics in torus winds. In a hypernova, these wind energies deposit their 
momenta onto surrounding remnant envelope, providing enhanced kinetic energies to matter 
ejecta. Wind energies determined from observational data on kinetic energies in GRB- 
supernova events (e.g., (Reeves et al. 2002)) and their remnants hereby provide a means 
for constraining the gravitational wave-frequency, and circumventing uncertainties in model 
parameters. 

The energetics of black hole-torus systems in suspended accretion depends on the mass 
of the black hole and the angular velocity of the torus. Hypernovae are believed to involve 
stellar mass black holes of about A — IAMq, consistent with the observed mass distribution in 
soft X-ray transients. The angular velocity of defines the efficiency of the catalytic conversion 
of black hole spin-energy as well as the frequency of gravitational radiation (at various 
multipolc moments). Yet, it is difficult to determine from first principles, and it becomes 
of interest to consider ways to constrain the angular frequency of the torus from the output 
in torus winds. The specific approach presented here considers the the phenomenology of 
torus winds, as well as a detailed calculation on the energy output in gravitational radiation, 
winds and thermal and neutrino emissions. 

The fate of the torus wind energies is determined by its nearby environment. We here 
propose that the torus wind energy is largely deposited within and onto the surrounding 
remnant stellar envelope, in the form of radiation and kinetic energy, respectively. We thus 
attribute the X-ray line-emission energies in matter ejecta in GRB-supernova associated 
events (e.g., (Reeves et al. 2002; Ghisellini et al. 2002)) to irradiation from within, as the 
expanding envelope reaches optical depth of order unity. We note that the deposition of torus 
wind energies is essentially that of a point explosion, given the relatively short durations of 
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activity at hand (about 20 s in association with a long gamma-ray burst). 

It has further been suggested that hypernovae may be associated with supershells (Efre- 
mov et al. 1998; Kim et al. 1999) which are X-ray bright (Wang 1999). Indeed, these X-ray 
bright features appear to be distinct from collective supernovae (Bunne et al. 2001). A bary- 
onic component to torus winds represents ejecta, which may further account for chemical 
enhancements in the companion stars in some of the soft X-ray transients (Brown et al. 
2000), e.g., GRO J1655-40 (Israehan et al. 1999) and V4148Sgr (Orosz et al. 2002). 

A GRB-hypernova association is receiving increasing support from the association of 
long gamma-ray bursts with supernovae, as in notably in GRB 980425/ 1998bw (Galama 
et al. 1998) and, more recently, GRB 011121 (Bloom et al. 2002), GRB 011211 (Reeves 
et al. 2002) and GRB 020405 (Price et al. 2002). Indeed, long GRBs are associated with 
star- forming regions (Bloom et al. 2002), and hence possibly with young massive stars in 
binaries. This appears to be consistent with a host environment in molecular clouds (Price 
et al. 2002). In turn, this suggests closer consideration of an association of hypernovae and 
their remnants in molecular clouds. 

Calorimetry on systems consisting of a torus in suspended accretion around a rapidly 
spinning black hole therefore can be approached as follows. First, the beamed output in true 
GRB energies Ej has been determined to be a few times 10^°ergs (Frail et al. 2001), which 
is indicative of ultrarelativistic baryon poor jets of somewhat larger energies. Second, the 
energy Ey^ in torus winds can be determined from energies associated with matter ejecta and 
their remnants in the host environment, the impact on a binary companion star in a remnant 
soft X-ray transient, and from beaming in GRBs when derived from collimation by winds. 
Third, the energy emitted in gravitational radiation of Eg^ ~ 6 x lO^^erg for a ten solar mass 
extreme Kerr black hole may be determined by upcoming gravitational wave experiments. 
A fourth energy deposition in MeV neutrinos is ulikely to be detectable in the foreseeable 
future, given the relatively large distances of these transient events. The energies E^ and 
Egu, are both produced by the torus and are a function of its angular velocity. This suggests 
determining in an effort to predict Egy^ and hence the frequency fg^ of gravitational 
radiation, as Eg^; and fg^ are closely related. 

We shall derive simphfied expressions for the fractions Egy,/Erot-, Ey,/Erot and in dis- 
sipation Ediss/Erot, relative to the rotational energy Erot of the black hole. These results 
are used to define an observational constraint on the torus-to-black hole angular velocity by 
the energy emitted in torus winds. This serves to circumvent otherwise uncertain physical 

parameters, notably the viscosity in the torus which provides coupling between the inner and 
the outer face. We propose to determine E-^j from calorimetry on hypernova ejecta and their 
remnants. Accurate prediction of the gravitational wave spectrum from wind energies will 
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be useful in designing dedicated experiments for searches for gravitational radiation from 
hypernovae. 



2. Multipole moments in wide tori 



Quite generally, a torus tends to develop instabilities in response to shear. This can 
be studied analytically in the approximation of incompressible fluid about an unperturbed 
angular velocity 



M^l'^ja^l'^ (Papaloizou & Pringle 1984). 



where q > 3/2 denotes the rotation index and ^It — M^^'^/a'^^'^ 
In the inviscid limit, we are further at liberty to consider the evolution of irrotational modes 
in response to initially irrotational perturbations to the underlying flow (vortical ii q ^ 2) 
by Kelvin's theorem. This approach shows the Papaloizou-Pringle instability (Papaloizou & 
Pringle 1984) to also operate in wide tori (van Putten 2002). The neutral stability curves 
of the resulting buckling modes can be described in terms of the critical rotation index 
Qc — Qc{b/ci,m) as a function of the minor-to-najor radius b/a. Quadratic fits to these 
stability curves are 



qc{b/a,m) = < 



2A9{b/ay + 1.73 


(m = 


2) 


6A7{b/ay + 1.73 


(m = 


3) 


< UAib/af + 1.73 


(m = 


4) 


20.0{b/ay + 1.73 


(m = 


5) 


^ 0.85m2(6/a)2 + 1.73 


(m > 


5) 



(2) 



Instability sets in above these curves, stability below. We note that for m = 2 the critical 
value qc — 2 obtains for 



b/a = 0.3225, 



(3) 



associated with the Rayleigh stability criterion for the azimuthally symmetric wave mode 
m — 0. For large m, we use the numerical result of critical values b/a — 0.56/m for q = 2. 

A quadrupole buckling mode radiates gravitational waves at close to twice the angular 
frequency of torus (van Putten 2002). Because the buckling mode represents an internal flow 
of energy and angular momentum from the inner to the outer face of the torus, in which 
total energy and angular momentum of the wave remain zero, it is not subject directly to the 
Chandrasekhar-Priedman-Schutz instability. Nevertheless, it is stimulated by gravitational 
radiation-backreaction forces, as derived in the approximation of the Burke-Thorne potential 
(van Putten 2002). 
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3. Energy output from a torus in suspended accretion 



The suspended accretion state is described by balance of energy and angular momentum, 
pertaining to gravitational radiation, Poynting flux-dominated winds, thermal and neutrino 
emissions (van Putten 2001): 



where denotes dissipation in thermal and neutrino emissions, Jl^ = (^+ + fi-)/2 is the 
angular velocity of the torus, defined as the mean of the angular velocities of the inner 
and and outer faces with angular velocities Q±. The inner and outer faces are subject to 
Maxwell stresses, giving rise to torques t+ (ingoing) and r_ (outgoing). In what follows 2TrA 
shall denote the active magnetic flux supported by the torus, representing open magnetic 
field-lines connected to infinity and the horizon of the black hole. These are endowed with 
no-slip/slip boundary conditions, and thereby assume the angular velocity of the torus. The 
remainder of magnetic field-lines are inactive, making up a toroidal 'bag' of closed field-lines, 
both on the inner and the outer face, similar to the closed field-lines in pulsar magnetospheres 
(Goldreich & Julian 1969; van Putten 2001). Thus, t_ = A^f^Q^ associated with a fraction 
/k, of flux in torus winds to inflnity, and r+ = {^h — ^+)^^/ff associated with a fraction fn 
of torus winds entering the black hole. With this definition, fw + fn — 

We consider shear stresses and the resulting dissipation in the torus to be due to tur- 
bulent magnetohydrodynamical fiow, in response to competing torques acting on the inner 
and outer face. These competing torques promote super-Keplerian and sub-Keplerian mo- 
tions, and hence the torus assumes a geometrically thick shape. By dimensional analysis, 
we consider the constitutive relation = ^Af[Q,+ — Q_)^ with 7 a factor of order unity 
and Ar = Rh < Bf. >^/^ denoting the root mean square of the radial fiux averaged over the 
interface between the two faces with contact area 2TiRh. The detailed structure of the mag- 
netohydrodynamical fiow which develops determines the ratio '-)A^/A^, which parametrizes 
the effective viscosity. The constitutive law for P^ thereby contains one free parameter. In 
what follows, we parametrize the relative strength of the net radial magnetic field as 



where a and h denote the major and the minor radius of the torus, respectively. This choice 
of parametrization (5) is such that z becomes independent of the aspect ratio 6/a, whenever 
the magnetohydrodynamical fiow develops a fiat infrared spectrum up to the first geometrical 
break m* = [a/b] in the azimuthal wave- number. 




(4) 




(5) 
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A leading order expansion in the minor radius b of (1) obtains Q± = f^r(l ± S) with 
fla = denoting the mean angular velocity of the torus and Jl± the average angular velocity 
of the inner and the outer face. Here, 

denotes the slenderness ratio of the torus. In particular, we have [Q] ~ Q+ — Q_ ~ q^xb/a. 
The equations (4) can be now be solved for 

and the energy output Eg^^ in gravitational radiation. Explicit expressions for the energy 
output E^i^, and Ed in winds and the combined dissipation in heat and neutrino emissions. 
The model parameters are, therefore, the rotation index q, the dimensionless magnetohydro- 
dynamical viscosity z, the slenderness ratio 5, the ratio of angular velocities rj, in addition 
to flux fractions fn and fw 



3.1. Estimate of rj 

The first equation of (4) may be used to eliminate r_|_ in the second equation, flTrad 
fl+T+ — f2_r_ — Pd, which obtains 



(8) 



With the constitutive ansatz for Pd given above, it follows that Trad = 2A^[r2] — 2r_. The 
luminosity in gravitational becomes ^ QT^rad, in view of the fact that the lower order 
multipole moments are essentially in corotation with the torus. Here, 



rad 



A^ 



M 



-2/, 



(9) 



where a — 2qz — 2/^ by = qb/a. With forementioned expression for t+, the first 

equation of (4) with (9) obtains {rin — ^+)/ff = fw^- + Oi^- Writing D,± — Qt ^ we 
obtain 



77 = 



(10) 
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3.2. Estimate of Egyj 

An expression for the luminosity Lg^ in gravitational radiation may be obtained, upon 
noting that most of the black hole luminosity is irradiated into the torus, i.e.: Lh — f^+r+. 
(Output along an open flux-tube by the black hole in association with a GRB is subdominant; 
van Putten & Levinson (2002).) We write fl^ffj = 7]VLnHpHA^ = V {^t/^+) Vl+T+Qn/i^H - 
for substitution in (9) after multiplication by a/ fjj. By (10) with = f^r(l + 8), we have 
Lgy^/Ln = a-^flri/il + 5 - 77(1 + 5f), i.e.: Lg^/Ln ^ a/{a + + 5{a + + fl)), upon 
neglecting terms of the order r]5 and higher. The total rate of conversion of rotational en- 
ergy of the black hole, including dissipation in the horizon (Thorne et al. 1986), is given by 
Eh = flHi^H — ^+)/^^^ over a duration T = Erot/En- The ratio of energy released in 
gravitational radiation as a fraction of the rotational energy Erot, therefore, becomes 



4. Energy estimates for symmetric flux distribution 

Consider a symmetric flux-distribution, given by equal fractions of open magnetic flux 
on the inner and the outer face: 

fH^L = 1/2. (12) 

Writing (10) in the form of 77 = ffj/{2qz + (1 + 5)(/l — /^)), it reduces to f] = l/8qz. A 
fiducial value of 1] = 0.15 obtains by considering the Rayleigh value q = 2 and z = 0.56, 
which attributes the radial modes in the infrared magnetohydrodynamical spectrum to an 
equipartition between the unstable Papaloizou-Pringle buckling modes. The estimate (11) 
becomes 

Egnj _ 8qz - 2 

Erot ~ Sqz{8qz{l + S/4) - 1) 

The energy released in winds satisfies Eyj — a~^f^{l — SyEgy^, whereby 

E^ _ {\-5f 
Erot 8qz{8qz{l + 5 /A) - 1) 

Finally, (8) shows that the energy radiated by dissipation in the torus satisfies 

Ed _ ^Egu, ^ 25 Eyj 

Erot Erot 1—5 Erot 



(13) 



(14) 



(15) 
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The case of strong viscosity (large z) and small 5 reduce these expressions further to the 
three fractions 

Egyj/Erot ~ V: Ey,/Erot ^ Ej E^ot ^ ^5. (16) 

In association to GRBs, wc attribute a minor output of black hole-spin energy with 
the formation of a baryon poor jet of energy Ej along an open magnetic flux-tube of the 
black hole. A universal horizon half-opening angle defines Ej ~ O'^ to be standard, and 
a generally small fraction of rotational energy of the black hole (van Putten & Levinson 
2002). The simplest geometrical relationship obtains when 6h is given by the poloidal 
curvature of the magnetic field, i.e.. Oh ~ 10°(6M/a). Without fine-tuning, this gives a 
small fraction Ej/Erot ~ 10^^ released in baryon poor outfiows, accompanied by a spread 
of one order of magnitude in response to a spread in torus radius by a factor of two - in 
agreement with the observed GRB energies E^ (Frail et al. 2001). The associated wind 
energies Ey,/ E^ot oc {M/af vary somewhat less. 



5. Frequency estimate of quadrupole radiation 

The expressions (16) provide a link between the energy Eyj in torus winds and the 
frequency in quadrupolar gravitational radiation, set by the angular velocity 7] of the torus. 
We propose to estimate 7] from E,^, to circumvents uncertainties in qz. This gives a frequency 
of quadrupole gravitational radiation given by fgw — l/{2nM){Eyj/Erot), or 

The wind-energy scale of 4 x lO^^erg corresponds to 77 = 0.1 and M — 7Mq. We can 
refine (17) by including a factor (1 — 2r])/{l — 5/2), in view of (13) and (14). Higher order 
multipole moments are present for a torus of increasingly smaller width (see §2), which emit 
at commensurably higher frequencies, and probably at somewhat lower intensities. 



6. Calorimetry on ejecta and hypernova remnants 

Torus wind energies impact the remnant stellar stellar envelope from within. This 
results in enhanced kinetic energies in matter ejecta in GRB-supernova associated events 
and, ultimately, impact the host environment. For GRB 011211 (Reeves et al. 2002), we 
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have 

^ me,/3e, ^ 2 X 10^2 1^^^ erg (18) 

for ejecta of mass niej with half-opening angle 9ej and (3ej — Vej/c, where Vej — 0.1c is the 
velocity of the ejecta in terms of the velocity of light c. At present, the angle 9ej is not 
constrained observationally. It may be larger than the GRB beaming angle, in which case 
we should see an over abundance of hypernovae with dim or no GRBs. We remark that 
only beamed outflows - baryon poor from the black hole - have sufficiently high luminosiy 
density per steradian to plow through a stellar envelope (MacFadyen & Woosley 1998). 

As the envelope expands, its optical depth reaches order unity and is irradiated from 
within by the accumulated remnant radiation from the torus wind energies. We attribute the 
observed X-ray line-emissions to this continuum emission, whose energies are on the order 
of 10^2 ej.g (Ghisellini et al. 2002). 

Follow-up calorimetry may be pursued on hypcrnova remnants, possibly on X-ray bright 
supcrshclls. This could be pursued by comparing this population in a nearby galaxy with 
with the true GRB-hypernova event rate of about 5 x 10^^ times the supernova rate (Woods 
& Loeb (1998), including beaming), and identification of X-ray point sources (see Brown et 
al. (2000)). 



7. Conclusions 

We quantify the energetics of a torus in suspended accretion against a magnetic wall 
around a rapidly rotating black hole formed in hypernovae. The torus catalyzes most of the 
black hole-spin energy, radiating a major fraction in a band limited burst of gravitational 
radiation. The frequency sweep is about ten percent during the conversion of half the 
rotational energy of an extreme Kerr black hole. The durations should agree with the de- 
redshifted distribution of Tgo of long GRBs (Fig. 1 in van Putten (2002)). 

In the slender torus limit, simplified expressions are derived for energies emitted in grav- 
itational radiation, winds and dissipation in the torus. We attribute the narrow distribution 
in observed true gamma-ray energies to baryon poor jets produced in an open magnetic 
flux-tube on the black hole to geometry without flne-tuning: the fraction of rotational en- 
ergy released in baryon poor outflows is set by a horizon half-opening angle 9h given by the 
poloidal curvature of the magnetic wall. 

The frequency of gravitational radiation is strongly correlated to the energy output in 
torus winds by (17). We attribute the energies associated in matter ejecta in GRB-supernova 
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events to the deposition of torus wind energy within and onto the surrounding remnant stellar 
envelope. We propose to hereby constrain the expected frequency of gravitational radiation 
by calorimetry on these matter ejecta and their remnants in the host environment. The latter 
may be pursued by statistics on the kinetic energies in matter cjccta, and by calorimetry on 
hypernova remnants in association with molecular clouds, conceivably in the form of X-ray 
bright supershells. 
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